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Abstract. In this paper we propose to apply the analogy between gravity and 
condensed matter physics to relativistic Bose-Einstein condensates (RBECs), i.e. 
condensates composed by relativistic constituents. While such systems are not yet 
a subject of experimental realization, they do provide us with a very rich analogue 
model of gravity, characterized by several novel features with respect to their non- 
relativistic counterpart. Relativistic condensates exhibit two (rather than one) quasi- 
particle excitations, a massless and a massive one, the latter disappearing in the 
non-relativistic limit. We show that the metric associated with the massless mode 
is a generalization of the usual acoustic geometry allowing also for non-conformally 
flat spatial sections. This is relevant, as it implies that these systems can allow the 
simulation of a wider variety of geometries. Finally, while in non-RBECs the transition 
is from Lorentzian to Galilean relativity, these systems represent an emergent gravity 
toy model where Lorentz symmetry is present (albeit with different limit speeds) at 
both low and high energies. Hence they could be used as a test field for better 
understanding the phenomenological implications of such a milder form of Lorentz 
violation at intermediate energies. 
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1. Introduction 

Analogies in physics have often provided deep insight and inspiration to deal with 
fundamental problems. In the last 20 years, analogue models of gravity [1] had this 
role with respect to pressing issues in gravitation theory, such as the mechanism of 
Hawking radiation, the fate of Lorentz invariance at ultra-short distances and the nature 
of spacetime and gravity as possibly emergent phenomena. 

The general idea behind analogue models of gravity is that in many physical 
systems it is possible to identify suitable excitations that propagate as fields on a 
curved spacetime. In particular, there is a wide class of condensed matter systems that 
admit a hydro dynamical description within which it is possible to show that acoustic 
disturbances propagate on an effective geometry, the so called 'acoustic metric' [2]. In 
fact, the propagation of excitations in the hydrodynamical regime can be described 
through a relativistic equation of motion in a curved spacetime. So even starting from 
non-relativistic equations one can show that excitations are endowed with a Lorentz 
invariant dynamics where the Lorentz group is generally associated with an invariant 
speed coinciding with the speed of sound. Furthermore, these systems generically show 
a breakdown of such acoustic regime leading to Lorentz violations at high energies (as 
expected from the Newtonian nature of the fundamental equations). The richness and 
power of this analogy have led in recent years to several applications. 

There is nowadays a vibrant activity toward the realization in laboratories of 
analogue black holes [3] and toward the observation of the associated Hawking 
radiation [4]. Acoustic metrics showing black hole horizons are found in many of the 
different systems proposed to realize this analogy: phonons in weakly interacting Bose- 
Einstein condensates [5], Fermi gases [6], superfluid Helium [7], slow light [8], non- 
linear electromagnetic waveguides [9] and ion rings [10]. In [11] a possible mapping to 
transform the acoustic metric into a metric conformally equivalent to the Schwarzschild 
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one has been pointed out. There are attempts to create also some analogue of rotating 
black hole geometries, with the main aim to predict super-radiance of the modes [12], an 
amplification of the modes due to the presence of an ergoregion [13]. In this direction, 
however, the analogy is limited by the fact that the spatial sections of acoustic geometries 
are necessarily conformally flat.f Expanding universes can be simulated with such 
systems, but, for similar reasons, only de Sitter spacetime, or in general, flat (k = 0) 
Friedmann-Robertson- Walker (FRW) (see e.g. [1, 15] and references therein) spacetimes 
can be cast in the acoustic form. This experiment-directed work had also motivated a 
wealth of theoretical studies regarding the robustness of Hawking radiation against 
ultraviolet (UV) Lorentz symmetry breaking [16] and on its possible signatures [17, 18]. 

While the above subjects have represented what we might call the main stream 
of research on analogue models, in more recent years there has been growing attention 
on their application as toy models for emergent gravity scenarios [19]. The latter are 
driven by the idea that gravity could be an intrinsically classical/large-scale phenomenon 
similar to a condensed matter state made of many fundamental constituents [20]. In this 
sense, gravity would not be a fundamental interaction but rather a large-scale/number 
effect, which emerges from a quite different dynamics of some elementary quantum 
objects. In this sense, many examples can be brought up, starting from the causal 
set proposal [21], passing to group field theory [22] or the recent quantum graphity 
models [23] and other approaches (see e.g. [24, 25]). 

Analogue models, and in particular Bose-Einstein condensates, have represented 
in this sense a very interesting test field as they are in many ways ideal systems 
where the mechanism of emergence of symmetries, spacetime and possibly gravitation- 
like dynamics, can be comprehensively studied and understood. In particular, much 
attention has been devoted in the recent past to Lorentz symmetry breaking and its 
possible role in leading to viable mechanisms of emergent dynamics (see e.g. the related 
discussion in [26]). 

However, except for the very recent work on the Abelian Higgs model [27], 
generically the analogue models considered so far are characterized by non-relativist ic 
fundamental equations. This would imply, in the emergent gravity perspective, a 
transplanckian world characterized by a preferred system of reference and a Newtonian 
absolute space and time. In this work we are going to consider a new analogue gravity 
system, namely a relativistic Bose-Einstein condensate (RBEC) for which the above 
issue is not present, and we shall show that also in this case an emergent spacetime 
structure can be found. From the point of view of emergent gravity analogues, this 
system can then be studied as an alternative to the standard case of a transition from 
Lorentzian to Galilean relativity in the UV, showing instead a transition from Lorentz to 
Lorentz symmetry encoded in the different invariant speeds at different energy regimes. 
Moreover, these systems are richer than their non-relativistic counterparts as they are 

| In [14] a partial mapping of the analogue metric into a Kerr black hole geometry has been done. 
Being the spatial part of the acoustic metric conformally flat, the mapping was possible just for the 
equatorial slice of the Kerr metric. 
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characterized by two kinds of propagating modes, a massless and a massive one (which 
disappears in the non-relativistic limit): we shall show that just the massless one is 
described by a generalized Klein-Gordon equation in a curved spacetime. Finally, we 
will show that the analogy with gravity can be applied to RBECs also and that the 
resulting metric allows to produce and investigate novel classes of metrics, as for example 
the k — — 1 FRW metric that we propose in the final part of the paper. Moreover the 
spatial slices of the resulting acoustic metric are not conformally flat, possibly permitting 
much wider mapping. 

The plan of the paper is the following. In section 2 the relativistic description of a 
BEC is given. The equation describing the perturbations in such a system is described 
in section 3, and their dispersion relation is studied in section 4 in different regimes. In 
section 5 the emergence of an acoustic metric is shown. As an application, in section 6 a 
possible mapping of the relativistic metric into a FRW geometry with k = —1 is proposed 
and commented. We finally conclude with some remarks in section 7, discussing what 
lesson can be learnt from this model. 

2. Relativistic BEC 

Bose-Einstein condensation, i.e. the macroscopic occupation of a single state, may 
occur both for relativistic and non-relativistic bosons: the main differences between their 
thermodynamical properties at finite temperature are due both to the different energy 
spectra and also to the presence, for relativistic bosons, of anti-bosons. These differences 
result in different conditions for the occurrence of Bose-Einstein condensation, which 
is possible, e.g., in two spatial dimensions for a homogeneous relativistic Bose gas, but 
not for its non-relativistic counterpart — and also, more importantly for our purposes, 
in the different structure of their excitation spectra. In this section we briefly recall the 
thermodynamic properties of a relativistic Bose gas [28]-[33], discussing how the non- 
relativistic limit is obtained and comparing with the results for non- RBECs [34, 35]. The 
study of the excitation spectrum for a (generally moving and inhomogeneous) condensate 
is presented in the following section 3. 

The Lagrangian density for an interacting relativistic scalar Bose field 4>(x,t) may 
be written as 

t = - v#t . V - + v(t, *)) fi$ - um> \), (i) 

where V(t, x) is an external potential depending both on time t and position x, m is 
the mass of the bosons and c is the velocity of light. U is an interaction term and 
the coupling constants Xi(t,x) can also depend on both time and position. U can be 
expanded as 

Wt0;A,) = ^ / ) 2 + ^ / ) 3 + ---, (2) 

where p = (jy(j). The usual two-particle A20 4 interaction corresponds to the first term 
(A2/2)/) 2 , while the second term represents the three-particle interaction and so on. 
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Since the Lagrangian (1) is invariant under the global U(l) symmetry, it is possible 
to define the associated Noether current and conserved charge, the latter corresponding 
to iV — N, where N (N) is the number of bosons (anti-bosons). 

Then, for an ideal homogeneous gas, when both U and V vanish, the relation 
between the chemical potential \x and the temperature T = is given by [28] 

AT-iV=^K(/i,/?)-n fe (/i,/3)], (3) 

k 

where n^p, (3) = 1/ {exp[/3(\Ek\ — p)} — 1} is the average number of bosons in the state 
of energy \Ek\, with 

E 2 k = h 2 k 2 c 2 + m 2 c\ (4) 

Similarly, n k (p, (3) = l/{exp[(3(\E k \ + p)] — 1} is the corresponding number of anti- 
bosons. 

Introducing the number density n = (TV — N)/Q (where Q is the volume of the 
system), one obtains in d dimensions for non- interacting bosons the following relation 
between the critical temperature T c and the charge density n [28, 32]: 



:-1 sinh ((3 c mc 2 

cosh (/3 c \Ek\) — cosh (/3 c mc 2 ) ' 

where C = 1/ (2 d ~ 1 7r d / 2 T(d/2)) is a numerical coefficient. 

From (5) one can readily derive the non-relativistic and the ultra-relativistic limits: 
the former is obtained when /csT c <C mc 2 . In this limit the contribution of anti-bosons 
to (3) can be neglected (so n w N/Q) and one obtains 

2rf / n \ 2 ' d 

kBT < = — [am ) (6) 

(( denotes the Riemann zeta function), which is the usual result for the critical 
temperature of a non-relativistic ideal Bose gas [34, 35]. In the ultra-relativistic limit, 
ksT c 3> mc 2 , one obtains 

d _ x h d c d - 2 T{d/2)(2n) d 
(W) = 4m 7 r^r(ci)C(^-l) n - (7) 
From now on we focus on the case d = 3, where the homogeneous relativistic ideal Bose 
gas undergoes Bose-Einstein condensation. 

At T T c , when the relativistic bosons condense, it is then possible to describe 
the dynamics of the condensate at the mean-field level by performing the substitution 
0—7-0: the order parameter <fi satisfies then the classical equation 

' +V(t,x))cj> + U'(p;\ i (t,x))(j ) = 0, (8) 



c 2 dt 2 T V K 2 

where p = and ' denotes the derivative with respect to p. The nonlinear Klein- 
Gordon equation (8) gives the dynamics of the relativistic condensates, and in then 
non-relativistic limit the Gross-Pitaevskii equation is retrieved. 
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Adopting the standard definition for the box operator in flat spacetime, 

□ = TTdfiv = -\dl + V 2 , (9) 



c 



equation (8) can be written as 



□0- [^ + v)<j>-u'{ P -\ i y 

3. Analysis of perturbations 

In this section, we study the excitation spectrum of perturbations on a condensate 
obeying the classical wave function equation (10). The field can be written as a 
classical field (the condensate) plus perturbation: 

= 0(1 + ^). (11) 

It is now worth noting that the expansion in (11) can be linked straightforwardly to the 
usual expansion [1, 34, 35] in phase and density perturbations 6\, p\\ 

hj_±^t * 1 = izil. (12) 

p 2 ' 2i v ; 

The equation for the quantum field ip describing the perturbations is 

U^ + 2r l ^{d^<P)d^- P U"(p-\ i )^ + ^) = 0. (13) 

It is now quite convenient to adopt a Madelung representation for the complex mean 
field and decompose it into two real fields, its modulus p(x, t) and its phase 6(x,t) 

= Vpe id , (14) 

the logarithm in (13) then becomes 

d ll \nct>= l -d ll \Yi P + id li e. (15) 
For convenience we define the following quantities: 

u m = —rj^duO, (16) 
m 

4 = ^pU"(r,>*), (17) 

T, = (□ + rTd^pdu) = -^f^4, (18) 

where the derivatives act on everything on their right, Cq encodes the strength of the 
interactions and has dimensions of a velocity and T p is a generalized kinetic operator 
that reduces, in the non-relativist ic limit c — > oo and for constant p, to the standard 
kinetic energy operator T: 

T P -)• -^-VpV = -^V 2 = T. (19) 



Imp 2m 
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A straightforward physical interpretation can be given to the four- vector u^. One can 
introduce the conserved current 

J, = 7^7 ( W - 0*^0) , (20) 
and show that 

Tfi 

j» = pdfl = P-^Uf,, (21) 

hence relating to the current associated with the U(l) symmetry. 

In terms of these quantities and using the phase-density decomposition, the 
equation for the condensate classical wave function, (10), becomes 

d^pu") = 0, (22) 



h 



,2 

-u = c 2 + 



V(x») + U'(p; Ai(a^)) 



(23) 



y/p 

The first equation is a continuity equation, which tells that the current j M defined above 
is conserved. The second one allows one to determine the zero-component of w M and, 
equivalently, the chemical potential, as a function of the spatial part of the fluid velocity, 
the strength of the interaction and the condensate density p. 

Multiplying (13) by h 2 /2m, one can rewrite the equation for perturbations in an 
RBEC as 

[ihu^dp -T p - mc 2 ] $ = mc ft. (24) 

Taking the Hermitian conjugate of this equation, we can eliminate ip\ obtaining a single 
equation for ip: 

[-ihu^ -T p - mc 2 ] [ihu^dp -T p - mc 2 ] $ = m 2 c\ ip, (25) 
and, with some simple manipulations, we obtain 

{[ihu^dp + T p ] 1 [-iku^dp + T p ] - -^dppdA $ = 0. (26) 
I c o P J 

This is the generalization to a relativistic condensate of the equation describing the 

propagation of the linearized perturbations on top of a non-relativistic BEC, in the 

same form as in [18]. It is worth stressing that (26) is implied by (24) but the converse 

is not true. Nevertheless (26) is enough for our aim, since we are interested only in 

the dispersion relation and in the metric describing the field propagation, which can be 

obtained from the modified Klein-Gordon equation in a curved background. 

As an application of (26), we consider a moving condensate with homogeneous 

density, V(t, x) — and A, being constant both in space and time. Let us choose 

0(t,x) =<j) e i{q - x ^ lt/h \ (27) 



one has 



H = mcu°, (28) 
q = mu/h, (29) 
= {n/mc 2 )d/dt + (h/m)q ■ V. (30) 
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Moreover (23) reduces to 

^ = h 2 c 2 q 2 + m 2 c 4 + h 2 c 2 U'(p ; A;), (31) 
where p = 4>o<Po- 

Note that p is the relativistic chemical potential of the condensate, which is related 
to the non-relativistic counterpart /inr via the identity /i = mc 2 + /inr- In particular, 
from (31), in the non-relativistic limit /inr <C mc 2 and fi ~ mc 2 . In the general case, 
when the fluid is not homogeneous, the same argument implies, from (23), that u° « c. 

Setting ip oc exp[i (fc ■ x — cut)] one obtains from (26) 

—u ■ k H uj -uj 2 H /c 2 

\ c 2mr 2m J 

X Lfe_i_^ _ 2 fc2 = a (32) 

V c 2mc 2 2m / V c / V ' 

For a condensate at rest (u — 0), the spectrum (32) was derived in [28, 33, 36]. 

The non-relativistic limit of (26) can be recovered by letting c — > oo. This implies 
/i = mc 2 + /inr ~ mc 2 and, equivalently w° ~ c. In this limit, the mode equation 
becomes 




Equation (33) reduces to the usual equation describing the propagation of perturbations 
in a non-relativistic BEC. 



4. The dispersion relation 



Before showing how the propagation of phonons in an RBEC can be described 
through an effective metric, it is worthwhile to analyze the dispersion relation of such 
perturbations. To do that we assume in this section that u, fi, p and Co are constant 
both in space and in time. For simplicity we start to analyze the case of background 
fluid at rest, u = 0. As discussed after (26), in order to study the dispersion relation, 
one can directly work on (26), which becomes in this case 



c 2m 



c 2m 



c 2 D 



^ = 0, 



(35) 



which can be rewritten expanding the □ operator as 

h w 2 ^ 



c 



2mc 2 



x 



2m 



-i—dt 

c 



h 



2mc 2 



Of 



w 2 

2m 



7)« 



clV 2 



i> = o 



(36) 
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The previous equation can be solved exactly by Fourier modes exp(— iut + ik ■ x) 



u 



-UJ 



c 2mc 2 
whose solution is: 

cui = c 



2m 



mu 



u 

c 2mc 2 



Co 



to 



(1)^ + ^ = 0, (37) 




mu y 



(35 



The above equation represents the dispersion relation for the modes in an RBEC 
and it is the generalization of the non-relativistic Bogoliubov dispersion relation. 

One can see from (38) that uf± ^ 0, i.e. there is no dynamical instability: in figure 1 
we plot the dispersion relation (38) for the interaction U(p) = J2 for different values 
of the dimensionless parameter A = (h/mc) 2 \2Po- In figure 1, we also plot the gap at 
k = between oj + (k = 0) and u)-{k = 0). 

One sees from (38) that the different regimes allowed for the excitations of the 
systems are determined by the relative strength of the the first two terms on the right- 
hand side of (38) (note that the same terms enter in the square root). Hence, in what 
follows we shall analyze the RBEC excitations dispersion relation (38) in some significant 
limits. The results are summarized at the end of this section in table 1. 



4-1. Low momentum limit 



We begin our study by looking at the low momentum limit of the dispersion relation 
(38) characterized as 



\k\ < 



mu 



h 



co 



mu 



h 



(1 + 6), 



(39) 



where, for later convenience, we have introduced the dimensionless parameter 



(40) 



Under the above assumption (39), the square root in (38) can be expanded as 



4 



k +2 



mu 



0\ 2 



b) ±2 



mu 



h 



0\ 2 



;i + 6)± 



=F 



k 4 



A(mu°/h) 2 (l + by 



.(41) 



4-1.1. Gapless excitations. Let us focus for the moment on the branch corresponding 
to the lower sign in the above dispersion relation. One obtains 



u 2 « c 2 



h 2 J (AO) 

.1 + 6 4(mu°/H) 2 (l + 6) 3 J ' 1 ' 

This is the dispersion relation of a massless quasi-particle and has the same form of the 

Bogoliubov dispersion relation in a non-relativistic BEC. Let us study now when the k 2 

term dominates over the k 4 one and vice versa and check if the conditions obtained for 

these two regimes are compatible with (39). 
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Figure 1. Excitation spectrum (38) for a condensate at rest (u = 0) with interaction 
U(p) = \2p 2 /2 and A = (h/mc) 2 \2po = 0.01, 1, 100. Solid line: gapless branch w_, 
dashed line: gapped branch uj + . u> in units of mc 2 /h, k in units of rac/h. Bottom 
right: plot of the the gap Au; = lo + — w_ at k = 0, which in these units is given by 
Aw(fc = 0) = V4 + 3A. 



Phononic (IR relativistic) regime: It is easy to see that the quartic term fc 4 can be 
neglected whenever 

|fc|<2^(l + fo). (43) 

h 

This condition is always compatible with (39). Depending on the value of b, they can 
be written in the compact form 

|Jfe|<^(l + &)min{l,2^}. (44) 
In this limit the dispersion relation reduces to the usual phonon dispersion relation 

ul = c\k\ (45) 
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where the speed of sound c s is defined as 
(cc /n°) 2 
1 + (c /u°) 2 1 



cl - } CC :' U Z, = A- (46) 



Newtonian (UV Galilean) regime: The opposite regime, in which the k 2 can be 
neglected with respect to k 4 , is defined by the following condition: 

+ (47) 

This regime is present only in the low-coupling case when b <C 1, since only in this 
case (47) and (39) are compatible. They can be summarized as 

mc ... mu° 

— « i*i « — • < 48 > 

Under these assumptions the dispersion relation becomes 

= jmoL = (mo! (49) 

2m M °/c 2/j/c 2 ' { ' 

which represents the Newtonian dispersion relation for massive particles with effective 
mass mu°/c = fi/c 2 . This is not surprising, since it is known from non-relativistic 
BEC that decreasing the wavelength of the perturbation, the atomic structure of the 
condensate emerges. 

In this case (an RBEC) we see that, even if the structure of the Lagrangian is fully 
relativistic, when the energy of the perturbations is not much higher than the chemical 
potential, the bosons are moving with non-relativistic speed, and the dispersion relation 
is Newtonian. Noticeably, the mass of the non-relativistic quasi-particles is not the 
mass of the bosons m, but it is an effective mass fi/c 2 , the rest mass corresponding to 
an energy equal to the chemical potential /i. This is not surprising since the chemical 
potential actually represents the energy needed to add a particle to the condensate. 

The non-relativistic limit: We discuss the non-relativistic limit of the gapless dispersion 
relation (41) starting from the behavior of the parameter b = (co/u ) 2 . As shown in 
section 3, u° — > c in this limit when the fluid flows with a velocity much smaller than 
the speed of light. Moreover, in the same limit, the interaction must be weak, namely 
c <C c, such that 

b « (c /c) 2 < I. (50) 

In this sense we can say that b measures the relativistic nature of the condensate, taking 
into account both the strength of the interaction and the velocity of the fluid, as it is 
possible to see from the definition of b (40), using (17), (28) and (31). Using (50) in (41) 
it is easy to obtain the Bogoliubov dispersion relation for a non-relativistic BEC (see, 
for example, [1]) 

k 4 



u 2 _ «c 2 



k 2 + 



(2mc /hy 



(51) 
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It is worth noting that c s (46) is, in general, different from the speed of sound Co in a 
non-relativistic BEC, but, as it is evident from the above equation, c s reduces correctly 
to Co in the non-relativistic limit. The functional form of the dispersion relation (42) is 
unchanged with respect to (51). Therefore, at this energy regime, there is no hint of the 
relativistic nature of the underlying BEC, without knowing the microphysics governing 
the elementary constituents. In other words, at this energy scale, a phononic observer 
cannot be aware of the relativistic nature of the condensate. 

Finally, let us also note that in this limit, condition (44), defining the phononic 
regime, reduces to the usual 

l*| « 2^, (52) 

where mco/h is the inverse of the healing length. When this condition is satisfied, one 
obtains the phononic dispersion relation 

uZ = elk 2 , (53) 

while, in the opposite regime, the dispersion relation describes Newtonian particles with 
mass m 

h 2 k 2 , . 

^ = W ™ 
as expected taking directly the non-relativistic limit of (49) and putting \i w mc 2 . 

4-1.2. Gapped excitations. A very different situation appears for the upper sign solution 
in (41): 



u) 2 ps c 2 



mu°\ 2 , „ 2 + 6, 2 k A 
l + b) + - -k 



h J v ^ ' " y ' 1 + b'" 4{mu /h) 2 (l + b) 3 

The only regime allowed is now the one in which the k 2 term dominates over the k 4 
term.§ In this limit (55) becomes 

2 4 
7T? C 

4 = + (56) 
where we have defined 

2 + b 
1 



C s,gap — 1 , L ^ ' (57) 



(l + b) 3 / 2 u° (l + 6) 3 / 2 /i 

; fTl = 2 : 7, 

2 + b c 2 + b c 2 

Equation (56) is the dispersion relation for a massive relativistic particle of effective 
mass m e fj and limit speed Cs igap . Note that, when (39) holds, and for b 1, the k 2 term 
is always dominated by the mass term. That is, this mode represents non-relativistic 
particles with mass m G g. When instead b 3> 1, the k 2 term can be of the same order 

§ In fact, the latter term would dominate if k 2 > 4(mn°/S) 2 (l + b) 2 {2 + b)> 8(mu°/h) 2 (l + b) 2 , but 
the above condition is clearly not satisfied when (39) holds. 
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of magnitude as the mass term, and the fully relativistic dispersion relation with mass 
gap becomes important. The k 2 term dominates when 



Let us now spend a few words about the physical meaning of this massive mode. 
Let us assume, for simplicity, that b can be neglected. In this case, from the above 
formulae, the gap is Ao; ~ 2mc 2 /h. This is an excitation of completely different nature 
with respect to previously discussed phonons. The mass gap of 2m indicates that the 
lowest possible excitation in this mode requires us to create a boson-anti-boson couple. 
When k 7^ the massive excitation propagates: this can be seen as the creation of a 
couple of particles, then its annihilation, while its energy is used to create another couple 
close to the former one and so on. Note that, apparently, this mode can make energy 
travel faster than light. However, one should consider that the dispersion relation (56) 
is valid only for sufficiently low k, satisfying (39). It is easy to check indeed that under 
this condition for k, the group velocity du/dk always remains smaller than c. 

The non-relativistic limit of the gapped dispersion relation (55) is trivial. The zero- 
order term in k 2 becomes now 4m 2 c 4 /h 2 , which diverges when taking the limit c — > oo. 
In fact the gapped branch disappears because the energy needed to excite this mode is 
much larger (indeed infinite in the limit) than the typical energy scales in non-relativistic 
configurations. In other words, this mode cannot be excited because one cannot create 
particle-anti-particle couples in the non-relativistic limit. 

4-2. High momentum 

The situation for high momentum is much simpler. When 



it is easily found that the dispersion relation (38) assumes the very simple form 



This is the standard relativistic dispersion relation of a massless field propagating at the 
speed of light c. This means that, when the energy of the perturbation is much larger 
than the chemical potential, the Newtonian particle regime of (49) is overcome. The 
dispersion relation assumes a novel form that was absent in the non-relativistic case. 
Therefore at very high momenta, the perturbations can probe the relativistic nature of 
the background condensate. 

5. The acoustic metric 




(59) 




(60) 



2 2; 2 

u± = c k . 



(61) 



The propagation of phonons can be described in the formalism of quantum field theory 
in a curved background only when the relativistic quantum potential T p can be neglected 
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Table 1. Dispersion relation of gapless and gapped modes in different energy regimes. 
cl = c 2 b/(l + b), c 2 gap = c 2 (2 + 6)/(l+6), m eS = 2( M /c 2 )(l+6) 3 / 2 /(2 + 6), b = (c /u ) 2 . 







Gapless Gapped 






b < 1 6> 1 


\k\ < m«°(l - 


|fc|<2mc°/n 
h 0J 7 " 2mc°/h < | fe < 77iu° /Jl 


, ,2 _ _2 t 2 

— s , i 2 — ^ 2 t 2 /,, 2 — m 2 r- 4 /fi 2 4_r 2 h 2 
Hu, = (Hk) 2 /2{/J,/c 2 ) U ~ C ° k ~ m cff C s,gap/ ft + C S ,gap« 


\k\ > mu°(l - 


h 6)/fi 


^ 2 = C 2 fc 2 



in (26). || This is the relativistic generalization of neglecting the quantum pressure in 
the Gross-Pit aevskii equation in a non-relativistic BEC. 

This is possible under two conditions. Firstly, the frequency and wave-number on 
the gapless branch oj_, (42), must be sufficiently low (see (44)). Secondly, one has to 
perform an eikonal approximation, i.e. one assumes that all the background quantities 
vary slowly in space and time on scales comparable with the wavelength and the period 
of the perturbation, respectively. This latter condition is equivalent to requiring that 

(62) 

and the corresponding relations for the variation in space. Under these assumptions, 
(26) simplifies to 



dtp 


< u, 


d t c 


< OJ, 


dtUf, 




P 




c 









4 p 



1> = o, 



(63) 



where, again, the differential operators act on everything on their right. In order to 
recover the acoustic metric, we need to cast the above equation in the form 



where is the metric density 



0. 



/ 



/IV 



-99 



(64) 



(65) 



To this aim, as in the non-relativistic case, one can use the continuity equation (22) 

rTdpj v = d^pu") = 0, (66) 
to commute with <9 M 



PV 



d v i) = 0, 



from which the metric density is easily read 



r= P_ 



-ci - ( U y 

—u°u 



-u°u 



Cnl — U Cg> U 



(67) 



(6* 



|| As discussed after (26), to read off the dispersion relation and the metric describing the field 
propagation it is enough to start from (26) instead of the full (24). 
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Finally, the metric describing the propagation of phonons in an RBEC is 



9/u> 



P 



Urr'U 



(69) 



y/l ~ U a U a /c 2 _ 

in coordinates = (x°,x) = (ct,x). 

The acoustic metric for perturbations in a relativistic, barotropic and irrotational 
fluid flow has been very recently derived in [37]. The above metric can be put in the 
same form as in [37] by a few variable redefinitions. To achieve this aim, let us define a 
four- velocity v^: 

(70) 



-w 



\u\ 



11 



^J—u a u a . 



With this definition, it is possible to generalize (46) when the spatial part of the four- 
vector is different from zero. To keep all the formulas in a covariant form, the scalar 
speed of sound c s must be defined as 

„2 C*cl/\\uf 



Using the above definitions, the metric (69) reduces to 



= p- 



(71) 



(72) 



This metric is manifestly conformal to that found in [37]. However, it is worth checking 
that they are exactly the same metric, comparing the conformal factors too. We have 
to show that pc/c s matches perfectly with the hlc 2 /c s (po +po) in [37]. All the quantities 
of [37] will be tilded to avoid confusion with ours. From (54) of [37] we can write 

,po +Po 



\u 



C- 



where C is a constant and, by definition, 



\u\ 



(73) 

V6. Note that from (46) of [37] 



n 



\u\\p, such that, apart from irrelevant constant factors, 



r,2 



n 



P 



C s (po+Po) C s 

and also the conformal factors coincide. 



(74) 



6. An application: k — — 1 FRW metric 

Mapping the metrics describing expanding universes is one of the most interesting 
applications of analogue models. However, only a mapping of the de Sitter spacetime 
and the k = FRW has been obtained so far [1, 15].^[ 

% Since every FRW metric is conformally flat, the natural attempt to reproduce such a spacetime within 
the framework of the standard Eulcrian fluid-gravitational analogy is to set v = and c s = constant. 
However, with this choice only stationary universes can be reproduced, because the acoustic conformal 
factor p cannot depend on time by the continuity equation. Of course, one may try to map a non-flat 
FRW metric onto an acoustic one by allowing for timc-dependent v, c s and p. Nevertheless, this task 
(even though not forbidden by any no-go theorem) would be much more difficult than the natural 
procedure obtained by using a relativistic fluid. 
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Here we show how the above acoustic metric can be mapped into the FWR metric 
with k = — 1. This is just an application, but it is important to stress how the category 
of metrics that one can map in this acoustic form has been immediately generalized. 
We start by rewriting (72) in spherical coordinates, assuming isotropy and v e = 
= 0: 



= p- 



( -1 + A^f/c 2 
—Av t v r /c 2 





— Av t v r /c 2 
1 + A(v r ) 2 /c 2 





o o \ 



r 2 

r 2 sin 2 (#)/ 



(75) 



where 



A 



Defining the new coordinate r and £ as 



t = V c 2 t 2 — r 2 , 



and choosing the following velocity profile: 
V< = cy/l + e, 

v r = ct 

one can write the inverse transformations as 



(76) 



(77) 



(7? 



t = r 



v 



T- 



(79) 



r£ = r — . 



To investigate the physical meaning of these coordinates, it is worth applying 
the above transformation also to the underlying Minkowski spacetime, seen by the 
condensate. The fluid velocity v can be transformed to 



(80) 



v * = v 6 = v * = 0. 

Therefore they represent coordinates comoving with the fluid. The flat Minkowski line 
element dsM is 

d^ 

i + e 

Going back to the acoustic metric, with the above choices the acoustic line element 
ds 2 becomes 



ds 



M 



-c 2 dt 2 + dr 2 + r 2 dtt 2 



-dr 2 + r 2 



+ £ 2 d^ 



ID 



ds 2 = p- 



-^dr 2 + r 2 



(82) 



One may notice that, when both c s and p depend only on r, the above expression 
represents the line element of a FRW spacetime with hyperbolic space sections: 

de 



ds' 



—p(r) dr~ + r p{r)- 



C.AT 



1 + ? 



+ e art 



(83) 
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However p cannot be chosen freely because it is related through the continuity 
equation (22) to the velocity = u M /||u||, which has already been fixed in (78). 
Nonetheless, it is easy to show that assuming p = po constant both in r and £, one 
can satisfy the continuity equation with an appropriate ||u||. Using (80) and (81) the 
continuity equation simplifies to 

<9 r (r 3 |M|) = 0, (84) 



from which 



\u\\=c^, (85) 



where E is an arbitrary function. 

It is therefore possible to define a new time coordinate 



dT = A /po^dr (86) 



and a function a{T) 



a(T) = rJpo—-, (87) 
V c,(t) 

where r(T) is implicitly defined in (86). In this way, the line element assumes the 
familiar k = — 1 FRW form 

d£_ 2 

At least from a mathematical standpoint, one has the freedom to choose an arbitrary 
form for a(T). This is equivalent to choosing c s as an arbitrary function of r. Using (86) 
and (87) 



ds 2 = -dT 2 + a{T) 2 ( -^- 2 + f W J . (88) 




r(T) = \l- f a{T>)dT>, (89) 
V Po Jo 

from which one can get T as a function of r given a generic a(T), and c s (r) is 

*M - ^yr (90) 

Fixing c s (t) means fixing U(po; Aj(x M )) in (17), and hence Aj, as given functions of x M . 
Indeed, from (71), co(x M ) depends both on c s (r) and ||m(x m )||, the latter given by the 
continuity equation, as discussed above. Finally, going back to the Euler equation (23), 
for constant p, one can determine also \^(x M ), because all the other functions ||tt|| and U 
have already been fixed. Using a two-body interaction U(p; A,,) = \2(x tl )p 2 /2 one finds 



V{T,t) 



p h 2 C 2 — C 2 s {t) \ 

c 2 + c 2 (r) /3(0 N 



m 2 c 2 



1 + 



c 2 — c 2 (t) \ r 3 



(92) 



Note that for a given FRW spacetime and hence given a(T) and c s (r), one still has the 
freedom to choose the most convenient H(£). 
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7. Conclusions 

In this paper, we have applied for the first time the analogy between gravity and 
condensed matter to an RBEC. Starting from the relativistic description of Bose- 
Einstein condensates, we have studied the propagation of excitations on top of it in 
a very generic framework, allowing for non-homogeneous density and nontrivial flow 
velocity. The full dispersion relation for the modes is given and analyzed in different 
regimes (see table 1): there are two branches, a gapless one and one with mass gap. The 
former one has three different regimes: for low momenta (wavelength larger than the 
healing length in the non-relativistic limit) the dispersion relation is linear u 2 = c 2 s k 2 . 
Then, for intermediate energies, lower than the chemical potential of the condensate, the 
dispersion relation becomes that of a Newtonian massive particle, just as for the non- 
RBEC case, when the wavelength is larger than or comparable with the healing length. 
Finally, for an RBEC, a third regime appears when the energy of the perturbation is 
larger than the chemical potential, and the dispersion relation becomes linear again 
(u 2 = c 2 k 2 ), but with velocity equal to the speed of light. The second branch has 
instead a mass gap, showing two different regimes. In this case, the dispersion relation 
describes, for low momenta, a relativistic massive particle with some effective mass and 
an effective limit speed apparently larger than the speed of light. However, this limit 
speed cannot be reached by such perturbations because at energy exceeding the chemical 
potential, the dispersion relation enters the second regime, in which the frequency is 
linearly proportional to the speed of light c. In the non-relativistic limit this branch is 
inaccessible because the energy needed to excite such a mode diverges in this limit. All 
these features are evident from the plot of the two branches in figure 1. 

For sufficiently low momenta (a generalization of the non-relativistic condition 
involving the wavelength of the perturbations and the healing length), it is possible 
to describe the propagation of the phononic mode (massless branch) through the 
gravitational analogy, and the corresponding acoustic metric is given. This metric allows 
for a wider range of applications than the usual one obtained in non-relativistic fluids. 
We propose as an example the mapping to the k = —1 FRW metric for expanding 
universes. We also saw how it can be used to also map spacetimes whose spatial section 
is not conformally flat. This could give the possibility to map even more complicated 
metrics (such as for example the Kerr metric) or to build novel structures not allowed 
in the non-relativistic case. 

This application has another important feature. Differently from all the previous 
systems where the analogy with gravity has been applied, the system investigated 
here is fully relativistic. Therefore, it does not suffer from the troubles related to 
the breaking of the Lorentz symmetry due to the emergence at sufficiently short 
wavelengths of the Newtonian structure of the underlying fluid and of its constituents. 
Hence, the RBEC is the first example of emergent Lorentzian spacetime from a 
Lorentzian background, showing therefore a Lorentz-to-Lorentz symmetry transition 
at high frequencies. Transitions from an IR Lorentz symmetry to a UV Galilean one 
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are known to be severely constrained (see e.g. [38]-[40]). Furthermore, in effective field 
theories they do tend to 'percolate' from the UV to the IR via renormalization group 
effects, by showing up in non-suppressed corrections to the low energy propagators of 
elementary particles [41, 42]. We do wonder if situations like the one just exemplified by 
this novel analogue model of gravity might show a less severe behavior and ameliorate 
this naturalness problem common to most of UV Lorentz breaking theories. We hope 
to investigate this issue in future works. 
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